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Abstract 

We present evidence for the decay P^K^. The results are obtained from a data sample 

of 227 X 10*5 r(4S') BB decays collected with the BABAR detector at the PEPTI asymmetric- 
energy B Factory at SLAC. Prom a maximum-likelihood fit giving a yield of 99 ± 19 events and 
efficiency estimated from simulation we make a preliminary measurement of the branching fraction 
B{B^ p^K^) = (5.1±1.0±1.2)xl0-6 where the first error is statistical and the second systematic. 
The hypothesis of zero signal in the mass region, GOOMel^ — 930Mel^, is excluded at the 6.1a 
level. Allowing a B^ — > /o(600)-ftrg contribution in the fit allows us to exclude the hypothesis of 
zero B° p^K^ at the 3.5cr level. 
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Figure 1: The two main amplitudes expected to contribute to the decay B — > p Kg are shown 
above. These are a colour suppressed tree diagram (a) and a gluonic penguin (b). 



1 Introduction 

In the Standard Model (SM), CP violation arises from a single phase in the three-generation 
Cabibbo-Kobayashi-Maskawa quark-mixing matrix [1]. Any measurement indicating additional 
sources of CP violation would be evidence for new physics. A number of penguin-dominated 
decays [2] (states such as (pK^, r]'K^, K+R-R^, /o(980)K° and p^R°) offer potential for making 
such an observation: they carry the same weak phase as the decay J/tp R^ [3], neglecting 

CKM-suppressed amplitudes, and therefore the Standard Model predicts their mixing-induced CP- 
violation parameter to be —r]fsm2(3 = —rjf x 0.74 it 0.05 [4]. Heavy non-SM particles may appear 
in additional penguin diagrams, potentially leading to new C!P -violating phases and CP-violation 
parameters measurably different from those predicted by the Standard Model. 

The two leading diagrams for the channel B^ — > p^R^ are shown in Figure 1. The penguin 
diagram is expected to dominate. We present evidence for the decay B^ — > p^Rg . We take the 
quasi-two-body (Q2B) approach, restricting ourselves to the region of the 7r'''7r" -R'^, Dalitz plot 
dominated by the p^ contribution and taking effects due to the interference between the p^ and the 
other resonances in the Dalitz plot as systematic uncertainties. 

The data we use in this analysis were recorded with the BABAR detector at the PEP-II asymmetric- 
energy e'^e~ storage ring at SLAC. The data sample consists of an integrated luminosity of 205 fb~^, 
corresponding to (227 ± 2) x 10*^ BB pairs, collected at the T(4S') resonance ( "on-resonance" ) and 
16 fb"^ collected about 40MeV below the T{4S) ("off- resonance"). In Ref. [9] we describe the 
silicon vertex tracker and drift chamber used for track and vertex reconstruction, and the Cerenkov 
detector (DIRC), the electromagnetic calorimeter (EMC), and the instrumented flux return (IFR) 
used for particle identification. 



2 The Candidate Selection 



We reconstruct P^K^ candidates {B^^^ in the following) from combinations of a decaying 

to 'K^TT~ and a decaying to vr+vr^. For the 7r"'"7r~ pair from the candidate, we use the 
combined information of the tracking system, EMC, and DIRC to remove tracks positively identified 
as electrons, kaons, or protons. In addition, we require at least one of the two tracks to have a 
signature in the IFR that is inconsistent with the muon hypothesis. The mass of the p^ candidate is 
restricted to the interval 0.6 < rn(7r+7r^) < 0.93GeV/c^. To reduce combinatorial background from 
low momentum pions, we require | cos 9^^+ \ < 0.95, where is the angle between the directions 
of the positive pion and the parent B^ in the /?° rest frame. The if^ candidate is required to have 
a mass within 13MeV/c^ of the nominal mass [10] and a decay vertex separated from the p^ 
decay vertex by at least three standard deviations. In addition, the cosine of the angle between the 
Kg flight direction and the vector between the p° and the decay vertices must be greater than 
0.995. 

Two kinematic variables are used to discriminate between signal and combinatorial background. 
The first is the difference /S.E between the measured centcr-of-mass (CM) energy of the B candidate 
and \/s/2, where ^/s is the CM beam energy. The second variable is the beam-energy-substituted 

mass mES = \J {s/2 + pj • y*bY /Ef — p^, where the B momentum and the four- momentum of 
the initial T(46') state {Ei, pi) are defined in the laboratory frame. We require 5.23 < mss < 
5.29GeV/c2 and \/^E\ < 0.15 GeV. 

Continuum e~^e~ ^ qq {q = u, d, s, c) events are the dominant background. To enhance discrim- 
ination between signal and continuum, we use a neural network (NN) to combine five variables: 
the cosine of the angle between the B^^ direction and the beam axis in the CM, the cosine of 
the angle between the thrust axis of the B^^ candidate and the beam axis, the sum of momenta 
transverse to the direction of flight of the B^^^, and the zeroth and second angular moments Lo,2 
of the energy flow about the B^^^ thrust axis. The moments are defined by Lj = p^ x | cos 9i\^ , 
where 6i is the angle with respect to the B^^ thrust axis of the track or neutral cluster i, and 
Pi is its momentum. The sum excludes the tracks that make up the B^^ candidate. The NN is 
trained with off-resonance data and simulated signal events. The signal efficiency determined from 
Monte Carlo (MC) simulation is 25%. MC simulation shows that 13.8% of the selected signal events 
are mis-reconstructed, mainly due to combinatorial background from low-momentum background 
tracks being used to form the p^ candidate in place of one of the pions it decayed to. In total, 
21000 on-resonance events pass all selection criteria. 

3 Background from other B Decays 

We use high statistics MC-simulated events to study the background from other B decays. The 
charmless decay modes are grouped into eight classes. The six B^ decay modes to the TT~^TT^Kg final 
state are of particular importance since they have signal-like AE and iues distributions and their 
decay amplitudes interfere with the p^Kg decay amplitude. Among these modes are fo{980)Kg, 
f2{1270)Kg, K*^'it'~ (including other kaon resonances decaying to K^ir'^), and non-resonant B^ — 
'K~^'K~K^ decays. The inclusive charmless B^ 'k~^it~ branching fraction (23ib3) x 10~^ together 
with the available exclusive measurements [4], are used to infer upper limits on the contributions of 
these decays. Selection efficiencies are obtained from MC and used with these branching fractions 
to estimate the expected background. The charmed decay B^ D~Tr'^{D~ Tr~Kg) contributes 
significantly to the selected data sample despite the veto requiring the invariant mass of both Kgir 



Background Mode 


EfRciency (%) 


Brancnmg l^r action (10 j 


TV T 1 f* I 1 1 T 7 J 

JN umber or bxpected bvents 


ii:o*(1680)-7r+ 


0.04 ±0.01 


14± 14 


2±2 


SO ^ /o(980)ii:0 


0.40 ± 0.02 


2±2 


2±2 


K2*(1430) + 7r- 


0.20 ±0.02 


14± 14 


6±6 


^0 ^ /2(1270)i^0 


1.68 ±0.04 


1.7 ± 1.7 


6±6 




1.11 ±0.10 


2.8 ± 1.7 


7±4 


5° ^ Ko*(1430)+7r- 


0.13 ±0.01 


14± 14 


2±2 


B^ r]'K^ 


6.29 ±0.10 


6.6 ±0.9 


94± 13 




0.64 ± 0.03 


7.1 ±7.1 


10± 10 


BO ^ D-7r+ 


0.34 ± 0.02 


42 ±6 


32 ±8 


Charmed B^ 






171 ± 86 


Charmed B+ 






106 ± 58 



Table 1: Table of the B background modes included in the Maximum Likelihood fit. The first 
six modes in this table all decay to the final state tt^tt^K^ . The number of expected events 
refers to the full fit region. The error on the number of expected events (used to calculate the 
systematic error) is taken from the measured branching fraction or is set to 100% when the value 
is not obtained from a direct branching fraction measurement. 



combinations to be more than 40 Mev/c^ from the D~ mass as quoted in [10]. As a result it is 
dealt with as an individual component in the fit. Two additional classes account for the remaining 
neutral and charged B backgrounds, with b ^ c decays being the dominant component in both 
cases. 



4 The Maximum-Likelihood Fit 

We use an unbinned extended-maximum-likelihood fit to extract the p^Kg event yield. In view of 
a future analysis of time-dependent CP asymmetry, the events in this sample are flavour-tagged as 
B^ or 5° with the method described in [11]. Four flavour tagging categories and an "untagged" 
category are defined each having a different expected purity for the signal. The likelihood function 
for the Nk candidates in flavour tagging category k is 



riB 



n { ^seu (1 - /U)^^-^^ + /^R^^-^^ + NcuVf,, + Y: NB,e,,^S,k \ (1) 



where Nj^ is the sum of the yields of all components, signal and backgrounds, tagged in category 
k, Ns is the number of p^Kg signal events in the sample, is the fraction of signal events tagged 
in category k, f^-^ is the fraction of mis-reconstructed signal events in tagging category k, Nc^k 
is the number of continuum background events that arc tagged in category /c, and Nsjej^k is the 
number of i?-background events of class j that are tagged in category k. The i?-background event 
yields are fixed in the default fit. The total likelihood C is the product of the likelihoods for each 
tagging category. 

The probability density functions (PDFs) Pf"^^, P^~^^, and Vff,, for correctly recon- 
structed signal, mis-reconstructed signal, continuum background and 5-background class j (see 



Table 1), respectively, are each the product of the PDFs of four discriminating variables. Each 
signal and background PDF is thus given by: Vk = P(mEs) -^(AE^) •Pfc(NN) -Pd cos(6l^+)|) where 
mES) ^E, NNand | cos{6^+)\ are the variables described in section 2. 

The mES) AE", NN, | cos (^7^+)! PDFs for signal and B background are taken from the simulation 
There are 15 free parameters in the default fit, including the yields of continuum events in each 
of the tag categories and the signal yield. Nine of the free parameters are used to describe the shape 
of the continuum background - third order polynomials to describe the NN shape and | cos(0^+)|, a 
second order polynomial to describe the AE distribution and for uies ^ function describing phase 
space with a single free parameter for its slope. 

5 Fit Results 

The maximum likelihood fit results in a signal yield of 99 it 19, where the signal is the combination 
of p^Kl and 5° /o(600)K° events. When selection efficiency (25%), fraction of 

decaying to (50%), the fraction of decaying to tt+tt" (69%) and the initial number of B^ 
5° pairs are taken into consideration this leads to: 

B{B^ p^K^) = (5.1 ± 1.0 ± 1.2) X 10"^ 

where the first error is statistical and the second systematic. The hypothesis of zero signal in the 
mass region, 600Mey — 930Mey, is excluded at the 6.1(7 level. Allowing as a free parameter 
the yield of a S° — > /o(600)iir^ contribution in the fit, discussed in Section 6, allows us to exclude 
the hypothesis of zero B^ — > p^K^ at the 3.5cr level. 

Figure 2 shows distributions of AE, m-Es, |cos^^+| and m^+^-, that arc enhanced in signal 
content by cuts on the signal-to-continuum likelihood ratios of the other discriminating variables. 

We performed a number of validation fits with different fit configurations to test the stability 
of the nominal fit; the results are consistent within errors. We made use of simulated experiments 
where events are generated from the PDFs for signal, continuum and B Backgrounds (toy Monte 
Carlo). The likelihood of our sample was found to be in good agreement with that obtained from 
toy Monte Carlo. 

Among the additional tests we have performed are: 

1. Adding B^ —>■ p^K^ simulated events to data to check that the fitted yield increases with the 
number of added events. 

2. Fitting for B^ /o(980)i^g and observing a yield consistent with our previous analysis [5]. 

3. Fits where background yields are allowed to vary: this is attempted with B^ v'^s' charged 
and neutral charmed B contributions and /2(1270)i^3. Background yields are consistent with 
the estimated values used in the default fit, and signal yield does not change by a significant 
amount. 

6 Systematic Uncertainty 

The contributions to the systematic error on the signal yield are summarized in Table 2. To 
estimate the errors due to the fit procedure, we perform fits on 1000 toy Monte Carlo samples 
with the proportions of signal, continuum and 5-background events measured from data. A bias 
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Figure 2: Projections of (clockwise from top left) AE, tues, \ cos9{tt~^)\, m^+^- enhanced in p^Kg 
signal through a cut on the ratio of signal to background likelihood using all discriminating variables 
except the one plotted. The red (upper) curve represents a projection of the maximum-likelihood 
fit result. The blue (lower) curve represents the contribution from continuum events, and the green 
line (middle) indicates the combined contributions from continuum events and B backgrounds. 



Table 2: Summary of systematic uncertainties . 
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B-background 
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Signal Model 


4% 


Q2B Approximation 


4% 




20% 


Total 





of (3 lb 1)% is observed in these fits and the sum in quadrature of the bias and its error is assigned 
as a systematic uncertainty. We also perform similar fit tests using fully simulated signal and 
5-background events added in the correct proportions to toy Monte Carlo continuum events. We 
use this technique to pick up any biases that may escape the conventional toy test: we observe no 
additional statistically significant bias and hence do not add an additional systematic uncertainty. 

For each class of B Background, the expected event yields are varied according to the uncer- 
tainties in the measured or estimated branching fractions, and the change in signal yield taken as 
a systematic. 

The uncertainties due to the extraction of the signal PDFs from simulation are obtained from 
a control sample of fully reconstructed —>■ D~(—>- K^'k~)t:~^ decays. We assume that differences 
in PDFs between data and Monte Carlo for this mode imply an equivalent data/Monte Carlo 
discrepancy in B^ p'^Kg and use them to estimate a 4% systematic error. 

The systematic error introduced through the use of the quasi-two-body approximation, neglect- 
ing the interference effects between the and the other resonances in the Dalitz plot, is estimated 
from simulation. We use a toy Monte Carlo technique, allowing all relative strong phases to take 
random values in each experiment. We take the RMS change in the signal yield as the systematic 
uncertainty. All B^ decays to ir^ir^K^ expected to provide events that pass selection cuts are 
included in the study, including /o(980), /2(1270), and the K*~^n~ and higher kaon states. In ad- 
dition, a non-resonant 5° — > Tr+Tr^ifg component is included in the study. The proportion of each 
contribution is estimated using known exclusive measurements and the inclusive B° — Tr'^Tr''Kg 
rate. 

An important systematic effect is the possible presence of B^ /o(600)i^°, where /o(600) is 
used to denote a broad scalar contribution of ill defined width that may lie beneath the /9°. To 
determine an upper limit on the size of such an effect we perform a fit with a B° ^ /o(600)if° 
contributing a yield allowed to float freely (since we do not include vrvr mass in this fit we make no 
assumption about mass distribution). We fit a yield of 31 ± 27 5° ^ /o(600)iv:° events. We back 
this up by performing fits where helicity is not used as a discriminating variable in the full helicity 
range and in the \cos{6^+)\ < 0.5 regions. From the fit to the full range and the efficiency in the 
low helicity range (estimated from simulation) we expect 22 events in our data sample and observe 
36± 12. 

We take the difference between the default fit and a fit assuming a /o (600) Kg contribution 
estimated from this second study as an estimate of this systematic effect. This leads to a 20% 
systematic uncertainty in the signal yield. 



7 Summary 



In summary, we have presented evidence for p^Ks. Prom a maximum-likelihood fit of 

the signal yield and efficiency estimated from Monte Carlo we measure the branching fraction 
B{B^ p^K^) = (5.1 ± 1.0 ± 1.2) X 10"^ where the first error is statistical and the second 
systematic. The hypothesis of zero signal in the p^ mass region, 600MeV-930MeV, is excluded at 
the 6.1(7 level. Floating the yield of a i?'' ^ /o (600) Kg contribution in the fit, discussed in Section 
6, allows us to exclude the hypothesis of zero — > pK^ at the 3.5(7 level. 
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